Can immobile magnetic vortex nucleate a vortex— antivortex pair? 
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It is shown that under the action of rotating magnetic field an immobile vortex, contrary to a 
general belief, can nucleate a vortex-antivortex pair and switch its polarity. Two different kinds 
of OOMMF micromagnetic modeling are used: (i) the original vortex is pinned by the highly- 
anisotropic easy-axes impurity at the disk center, (ii) the vortex is pinned by artificially fixing the 
magnetization inside the vortex core in the planar vortex distribution. In both types of simulations 
a dip creation with a consequent vortex-antivortex pair nucleation is observed. Polarity switching 
occurs in the former case only. Our analytical approach is based on the transformation to the 
rotating frame of reference, both in the real space and in the magnetization space, and on the 
observation that the physical reason for the dip creation is softening of the dipole magnon mode 
due to magnetic field rotation. 

PACS numbers: 75.10.Hk, 75.40.Mg, 05.45.-a, 72.25.Ba, 85.75.-d 



A magnetic vortex forms a ground state of submi- 
cron sized magnetic particles, which provides high den- 
sity storage and high speed magnetic RAM [1]. One bit 
of information corresponds to the upward or downward 
magnetization of the vortex core (vortex polarity). Excit- 
ing the vortex motion by a high-frequency magnetic fields 
or by a spin polarized currents, one can switch the vor- 
tex polarity on a picoseconds time scale. It is known that 
the switching process is mediated by a vortex-antivortex 
pair creation [2] . Typically to observe this switching phe- 
nomenon, one has to excite the low frequency gyroscop- 
ical mode, which causes the vortex gyromotion. To the 
best of our knowledge the vortex core switching was al- 
ways observed experimentally and by simulations only for 
a moving vortex. Moreover, there is a strong belief that 
the switching occurs 'whenever the velocity of vortex- 
core motion reaches its critical velocity' Vcri [3], which is 
determined only by the exchange constant A] for typical 
soft materials is about Vcri ^ 300 m/s [3, 4]. 

In the present Letter we predict the switching for 
the immobile vortex by a rotating magnetic field. The 
switching picture also involves the mechanism of a dip 
formation followed by a nucleation of vortex-antivortex 
pair. We found that the driving force for the dip forma- 
tion is the in-plane magnetization inhomogeneity, cre- 
ated by a vortex together with the rotating magnetic 
field; the vortex out-of-plane structure and the vortex 
velocity are not principle for this mechanism. These con- 
clusions are confirmed by the micromagnetic simulations 
and by a simple analytical picture. 

Recently, we have reported about the vortex core 
switching by the homogeneous rotating field B{t) = 
{B cos cot ^ B sin cot ^0) in the ten GHz range [5], which is 
much higher than the gyrofrequency and is in the fre- 
quency range of the higher azimuthal mode with \m\ = 1 
[6, 7, 8]. Opposite to the pumping with the low gyrofre- 
quency, which results in the visible gyroscopic motion 



of the vortex position, the high frequency field leads to 
small amplitude (few nm) oscillations of the vortex po- 
sition, see the trajectories on Fig. 3 of Ref. [5]. Typical 
velocities when switching occurs are lower than reported 
in Ref. [3]. There appears a question: is it necessary for 
vortex to move in order to switch its polarity? 

To achieve switching of the immobile vortex we have 
performed numerically two different kinds of modeling 
using OOMMF micromagnetic simulations [15] 

(i) We have pinned the vortex by the highly- 
anisotropic easy-axes impurity at the disk center, see 
Fig. 1. Since the magnetization is always hold perpendic- 
ular to the disk plane within the impurity, shifting of the 
vortex core from the impurity leads to the increasing of 
magnetostatic energy of the system. Therefore one can 
use such type of impurity for pinning of a vortex with 
out-of-plane core structure in conditions of a weak exter- 
nal influence. Initially, the vortex has a positive polarity 
p = By applying the ac fleld, which rotates against 
the vortex polarity, oup < [in a clockwise (CW) direc- 
tion in our case], one can easily resolve the formation of 
the dip near the pinned vortex, see Fig. 1(a). Under the 
action of the field the dip deepens, see Fig. 1(b). When 
the amplitude of the dip reaches its maximum value, a 
vortex-antivortex pair is created, see Fig. 1(c). The posi- 
tions of the vortices and the antivortex can be identified 
by the cross-section of isosurfaces rrix = and ruy = 
[9]. One can identify from Fig. 1(c) positions of new born 
vortex and antivortex as well as the position of the ini- 
tial vortex, which is (as it is seen from the figure) still at 
origin. The further dynamics is well known [2, 9, 10, 11]: 
the new born antivortex annihilates with pinned original 
vortex, and eventually only the new born vortex with the 
opposite polarity survives (see Fig. 1(d)). 

(ii) In the second kind of simulations we artifi- 
cially pinned the vortex core by fixing the magnetiza- 
tion inside the vortex core in the planar vortex dis- 
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FIG. 1: (Color online). The modeling of the vortex-antivortex pair formation for the pinned vortex on the impurity in the 
Py disk (300 nm diameter, 20 nm thickness, 9 nm is a diameter of impurity with the easy-axis anisotropy K — 17 MJ/m^). 
Snapshots describe the magnetization dynamics at different times under the influence of the rotating in-plane magnetic field 
(a; = 6 GHz, B — 20 mT). The lower row shows the in-plane magnetization dynamics near the disk center. The in-plane 
magnetization distribution is indicated by streamlines with arrows. Green and yellow solid lines corresponds to the isosurfaces 
rrix = and ruy = 0, respectively. The out-of-plane magnetization profiles along the line ^ crossing the impurity is shown 
above the snapshots (right pictures). Left pictures indicate the in-plane magnetization dynamics in the whole disk. 




FIG. 2: (Color online). The vortex-anivortex pair formation in a Py disk (150 nm diameter, 20 nm thickness) within the fixed 
core modeling (see the details in the text) under the CW field influence (a; = 8 GHz, B — 2{) mT). Notations are the same as 
on Fig. 1. 



tribution, where the normalized magnetization m = 

^^^/T^^Tnf cos (/); — sin ^; m^^ takes a form = 

and (j) = X i '7r/2 with (r, x) being the polar coordi- 
nates in the disk plane. By this artificially pinned vortex 
state configuration, the central vortex is a priori immo- 
bile, and moreover, without out-of-plane structure. The 
direction of the dip magnetization is determined only by 
the direction of the field rotation: if we apply a CCW 
(CW) field, there appears a dip with > (m^ < 0), 
see Fig. 2. This high frequency field excites also magnon 
modes. One can resolve the mode with m = 1 on 
Figs. 2(a) and (c), and the mode with m = 2 on Fig. 2(b). 



During the pumping, the amplitude of the dip increases, 
and finally, the vortex-antivortex pair is created, see 
Fig. 2(d). Note that the new-born antivortex can not 
annihilate with the initial vortex in the fixed core model: 
it either annihilates with the new born vortex and the 
scenario repeats again and again, or it stands near the 
original fixed vortex, while the new born vortex is in- 
volved by the field. 

The goal of our study is to reveal the mechanism of 
the dip formation. We follow the simplest model of fixed 
core planar vortex. Numerically, we found that the dip 
formation occurs in well defined range of the field pa- 
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FIG. 3: (Color online) Diagram of the frequency range, 
where the dip formation occurs for a fixed diameter disks 
(2R — 150 nm) under the influence of a fixed field amplitude 
{B — 20 mT) rotating field. Colors indicate the depth of the 
dip. The circle over black symbol indicates that the vortex- 
ant ivortex pair creation occurs. Light-green symbols indicate 
frequencies of the azimuthal eigenmode with m = 1. 
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FIG. 4: (Color online) Amplitude of m(t) at the minimum 
(dip center) for the disk(2i? = 150 nm, h — lb nm) under the 
influence of rotating field with a fixed amplitude B — 20 mT 
and different frequencies. Numbers in rectangle frames de- 
notes frequencies in GHz. The inset contains the time de- 
pendence of the dip depth — m(t) obtained from the effective 
Lagrangian (2) for cjo = 2, cji — 1, 6e = 0.1, ?7 = 0.1, ci — C2 — 
0.25, C3 = 0.1. Corresponding frequencies in dimensionless 
units are shown in oval frames. 



rameters {uo^B)^ see Fig. 3. When the frequency of the 
field is less than a critical value c^ci, the dip is not formed. 
The same happens for the high frequency field {uo > ujc2)' 
More accurate, out of the borders of the range {ujd^ ^c2), 
the depth of the dip rapidly decreases, what is indicated 
by colors on Fig. 3. Both critical values cOci and u;c2 
depend on the disk thickness, which counts in favor of 
nonlocal magnetostatic nature of this phenomenon. An 
example of the temporal dependence of the dip depth is 
plotted on Fig. 4 for a certain disk thickness (15nm), for 
which a;ci=8GHz and a;c2=12GHz. For co < cOci just as 
for CO > U0c2 the dip depth does not reach its minimal 
value rriz = —1, unlike the case uOci < oo < {jJc2^ when 
the vortex- ant ivortex pair is born. The depth of the dip 
usually makes a number of oscillations while it reaches 
the limit value. Only for frequencies, which are close to 
the eigenfrequency of the mode with |m| = 1, it changes 
monotonously. 

The magnetic energy of the system under consideration 
consists of three terms: = J d^r (Wex + ^ms + ^f), 
where IVex = (Vm)^ is the exchange energy 
density, ly^s = ^/ dV'(m(r)-V) {m{r') ■V)\r -r'\-^ 
is the density of magnetostatic interaction energy. And 
= —b^/Y^^m^ cos{(l) — cot) is the interaction with a 
magnetic field b = {b cos ut^b sin uot^O). Here and be- 
low normalized quantities are used: = £^/(47rM^), 
b = B /{AttMs)^ CO and t are measured in units Q and 1/Q 
respectively, where Q = Air^yMg and 7 is gyromagnetic 
ratio. In the no-driving case b = the magnetic energy 
(o of cylindrical nanoparticles is invariant under transi- 
tion into a rotated frame of reference: x = X ~ 4^{x) = 
0(x) ~ ^- This means that the quantity J = + Lz 
with Mz = /d^r(l — m^) being the magnetization along 



the cylindrical axis z and = — J d^r{l — mz)d^(j) being 
the z-component of the orbital momentum, is conserved. 
When the ac field is turned on, the total momentum J is 
not conserved anymore but its existence allows to obtain 
that in the rotating frame of reference x = X~^t^ r = 
i = (^(x) = </^(x = X + ^t) — cot the magnetic energy is 
time-independent and has the form 



^ = ^ 



where the last term represents the rotation energy. The 
dynamics of the system is governed by the Landau- 
Lifshitz-Gilbert equations which in the rotating frame 
of reference have the form 



drriz 



-ml S4) 



0, 



5^ 2^ 



0, 



where ^ = (1 + 7/^) / d^x{mz — l)dt4> — ^ is the La- 
grangian and 77 is a damping constant. 

To gain some insight how the interaction with the mag- 
netic field (which is static in the rotating frame) together 
with the rotation provides the dip creation we use the 
Ansatz 



rriz = ao{t)fo{r) + [ai{t) cosx + Q^-i(^) sinx]/i(r), 
^ = X+ J+/5oW^o(r) 



(1) 



+ [Pi {t) cos X + P-i{t) sin x] 91 {r) , 



where {fm{r)^ gm{r))e'^^^ are the magnon eigenfunctions 
with the eigenfrequencies co^ {m = 0, 1) on the planar 
vortex background (see e.g. [12] for general case) and 
a^, and /3m^ are time dependent coefficients. Inserting 
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the Ansatz (1) into the Lagrangian ^ and carrying out 
integrations over the spatial coordinates, we get an effec- 
tive Lagrangian in the form 



m=0,±l 



\m\ 



m=0,±l 



+ cg/^o (^1 + + ^higher order terms^ . 



(2) 



Here the coefficients Cj with j = 0, . . . , 8 are due to non- 
hnear terms in the magnetic energy S' and be is an effec- 
tive strength of the magnetic field. Equations of motion 
for the parameters (5^ have the form 



V 



(3rn = mUJp-rn + ^ ^ 

dam 



dam ' 



(3) 



Thus the problem under consideration is reduced to 
the set of three nonlinear ly coupled oscillators. When 
uj <^ uji the linear terms in Eq. (2) dominate and 
the out-of-plane component of the magnetization m = 

^QfQ + Qf^ + a^_i is small, see inset in Fig. 4. However, 
as a result of rotation the magnon frequencies in the ro- 
tating frame of reference are shifted, ujm = oJm—'^^^ and 
the dipole magnon mode (i.e. the azimuthal mode with 
m = 1) becomes soft, which can be identified from Fig. 3. 
Near the threshold where uji is small, the nonlinear terms 
in the Lagrangian ^ become crucial and the system goes 
to a new regime with the finite value of m, see the in- 
set in Fig. 4. The appearance of the finite fh we identify 
with the dip creation. The dip direction is determined by 
the product 6e^- when beUJ > 0(< 0) the out-of-plane 
component of magnetization is negative (positive). This 
result is in a full agreement with the results of full-scale 
numerical simulations. 

In conclusion, we showed that in the presence of immo- 
bile planar vortex the rotating magnetic field produces a 
vortex- ant ivort ex pair. This process is the most efficient 
in a finite frequency interval which is determined by the 
aspect ratio and material properties of the nanoparticle. 
The physical reason of the dip creation with a conse- 
quent vortex- ant ivort ex nucleation is softening the dipole 
magnon mode and magnetic field induced breaking cylin- 
drical symmetry in the rotating frame of reference. In 



some respects such softness of modes due to rotation is 
similar to the problem of instability of BEC under ro- 
tation [13], and to the ZeFdovich-Starobinsky effect for 
a rotating black hole [14]. Under the continuous pump- 
ing, the soft mode is excited, and the system goes to the 
nonlinear regime, which results in the dip formation. 
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